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Abstract 

The standard operation of a dc SQUID leads to oscillatory electric fields that emit electromag¬ 
netic radiation and can change the state of the measured sample. A stationary SQUID could be 
advantageous when back action on the measured sample has to be avoided. We study a super¬ 
conducting loop that encloses a magnetic flux, connected to a superconducting and to a normal 
electrode, when a fixed electric current between the electrodes flows through the loop. The consid¬ 
ered circuit does not contain Josephson junctions. We find that in a very broad range of parameters 
the current flow converges to a stationary regime. The potential difference between the electrodes 
depends on the magnetic flux, so that measuring this voltage would provide information on the 
enclosed flux. The influence of thermal noise was estimated. The sizes of the voltage and of the 
power dissipation could be appropriate for the design of a practical fluxmeter. We found narrow 
ranges of flux at which the voltage varies sharply with the flux. 

PACS numbers: 85.25.Dq, 74.78.Na, 03.65.Ta 
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I. INTRODUCTION 


Superconducting quantum interference devices (SQUIDs) are the most sensitive available 
fiuxmeters.- - SQUIDs are classified into rf SQUIDs, which consist of a loop interrupted by 
a single Josephson junction and are monitored at radio frequency, and dc SQUIDs, with two 
Josephson junctions. The dc SQUID can really be used as a “dc” fluxmeter by gradually 
increasing the current through it, until its flux dependent critical current is reached; however, 
in the more practical procedure, the current through the SQUID is larger than its critical 
value, giving rise to a flux dependent average voltage. 

More in detail, the voltage in the dc SQUID oscillates at a frequency that follows from the 
Josephson relation. As a consequence the SQUID emits electromagnetic radiation, which 
exerts a back action on the measured sample. In the case of nanoscopic samples, this back 
action can alter the state that the SQUID was intended to measure. Therefore, there are 
situations in which a stationary instrument can be advantageous over the standard SQUID, 
even if less precise. 

Stationary fiuxmeters do exist, such as the ballistic Hall magnetometer,- the normal- 
metal-insulator-superconductor (NIS) interferometer,- or the superconducting quantum in¬ 
terference proximity transistor/ Here we propose a simple circuit, which under appropriate 
circumstances could be the fluxmeter of choice. 


II. PROPOSAL 


The circuit that we propose is as follows. A loop of superconducting material encloses 
the magnetic flux <f> that we intend to measure. Short wires at the left and the right connect 
the loop to “banks” that serve as electrodes. The connectors and the left bank are made of 
the same material as the loop, whereas the right bank is a normal metal. A fixed current is 
driven from the superconducting to the normal bank, and the potential difference between 
them is measured. If this potential difference is a function of <f>, then this circuit provides a 
measurement of $. Our proposal is sketched in Fig. [0 

Additional ideas for SQUIDs without Josephson junctions have been considered in 


Refs. 
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The difference between the circuit in Fig. |T]and the NIS interferometer is that 


the insulating barrier has been replaced with the connector to the normal bank. Another 
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FIG. 1: Circuit designed to measure the flux <I>. The circuit consists of a superconducting loop 
that encloses the flux, with connections to two electrodes. The electrode at the right is made of 
a normal metal. The cross section of the connector at the left is ws and the cross section at the 
right is u>n- Both branches of the loop have the same length. The cross section in each branch is 
the same function of the arc length, and varies monotonically from ws to u;n- 


circuit that can operate as a SQUID and involves normal components was investigated in 
Ref. [n|, but the the authors did not address the question of whether a stationary regime is 
attained. 

Heuristic arguments Due to fluxoid quantization, current will not only flow from one con¬ 
nector to the other, but a circulating current will also be present. It follows than the current 
in one of the branches is expected to exceed half of the total current, increasing the amount 
of current that has to be transported as normal current. We may therefore expect an in¬ 
crease in the voltage between the banks as the enclosed flux deviates from an integer number 
of quantum fluxes. This is true also for a standard SQUID above its critical current, but 
in the standard case the phase of the order parameter changes at a different rate in each of 
the banks, forcing an oscillatory behavior. In our case, the order parameter vanishes at the 
normal bank, so that the phase becomes meaningless, and oscillations will not necessarily 
occur. 

The current between the banks should not be too large, since in this case a considerable 
segment of the loop becomes essentially normal, making the circuit insensitive to the flux. 
On the other hand, the current ought to be sufficiently large to cause the order parameter 
to vanish at a point in the loop when the flux becomes half a quantum. This feature enables 
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a continuous passage between consecutive winding numbers, thus avoiding hysteresis. 


III. MINIMAL MODEL 


A. Formalism 


For simplicity, we consider a quasi-ID circuit, and the position dependence of all physical 
quantities will be completely determined by the arc length. In this section we assume that 
thermal contact is sufficiently effective to take heating away, so that the entire circuit is kept 
at a uniform temperature T; the expected influence of Joule heating and of supercurrent 
heat transport will be examined in Sec. IIVBI At this stage we will also ignore thermal 
fluctuations. 

In Refs. and 0 it is assumed that the entire voltage drop occurs across the insulating 
barrier, so that the superconducting circuit has a uniform potential and a static description 
is possible. In contrast, in our case the electrochemical potential varies along the super¬ 
conducting circuit, so that a time dependent formalism is necessary, and convergence to a 
stationary regime cannot be assumed a priori. 

This section will rely on time-dependent Ginzburg-Landau model (TDGL), which is the 
simplest model for the description of the proposed circuit. More realistic models will be 
considered in Sec. m 

The proximity effect is built-in into the Ginzburg-Landau model: near the superconduct¬ 
ing electrode, superconductivity will be strong, even above the critical current, and near the 
normal electrode superconductivity will be weak, even for very small currents. As a conse¬ 
quence, part of the current along the circuit will be normal current, giving rise to a potential 
difference between the electrodes. 

We choose a gauge with no tangential vector potential along the connectors and uniform 
tangential vector potential A along the loop. We write e = 1 — T/T c , where T c is the critical 
temperature of the superconductor, and denote by <p the electrochemical potential. The unit 
of length will be denoted by xq , by to the unit of time, by <po the unit of voltage, by A 0 the 
unit of vector potential, and by jo the unit of the current density. We take 
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Here £(0) is the coherence length at T — 0, ks is Boltzmann’s constant, e is the electron 
charge, and a is the normal conductivity. 

Let us first consider the case of uniform cross section, ta s = w N ; in this case the current 
density is the same at both connectors (assuming electroneutrality), and will be denoted by 
j. With the units in Eq. (IT]) , the TDGL equation and Ohm’s law read 

dti> = £> 2 0 + (e - |0| 2 - i<p) 0 (2) 


and 

d x tp = u 1m (0*1)0) — j . (3) 

Here 0 is the order parameter, with normalization imposed by Eq. (I2]h dt and d x denote 
partial differentiation with respect to the time t and to the arc length x along the wire, T> is 
the operator T> = d x —i A and u is the ratio between the relaxation times of 0 and j Along 
the branches of the loop, j has to be set as the current density that enters the corresponding 
branch. Following Ref. 12, we denote by 4L the perimeter of the loop, and therefore the flux 
that it encloses is 4LA Invoking the units in Eq. ([T]) we obtain |A| = 7r|<I>|/2L < f ) o, where 
4>o = 7T hc/e is the quantum of flux. Accordingly, along one of the branches of the loop we set 
A —>■ 7T<f ) /2L < f>o, along the other branch we set A —>■ — tiQ/2LQ 0 (for our purpose it makes no 
difference which branch is assigned the positive sign), and set A —» 0 along the connectors. 

The boundary conditions require continuity of the order parameter and assume equilib¬ 
rium at the electrodes. The potential at the superconducting electrode is taken as zero, and 
the order parameter as 0s = e 1 / 2 . At the normal electrode, the order parameter is required 
to vanish. At the contacts between the ring and the connectors, we require continuity of the 
potential and of the order parameter. We also impose the constraint that the sum of total 
currents along the branches of the loop has to equal the current along the connectors. 

Equations (|2 |) -(|3| ) are invariant under the gauge transformation A —y A + C, 0 —> 
exp (iCx)ijj where C is constant. If C is an integer multiple of n/2L, also the continuity 
conditions remain invariant. It follows that any physical property of the circuit is periodic 
in <3>, with periodicity $ 0 . In addition, switching the sign of $ amounts to exchanging 
the lower and upper branches of the circuit, so that the potential difference between the 
electrodes has to be a symmetric function of <f>. 

TDGL is not expected to be quantitatively correct, but one of the features that makes 
it a valuable tool is the scaling with length. Equations (l2l)-(|3|) are invariant under the 
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transformation x —» L'x , L —> L'L, A —>- L'~ x A , t —> L /2 t, tj} —>■ Z/ -1 ^, </? —> L'~ 2 ip, 

e —)■ L /_2 e and j —>■ L ,_3 j, since each of the terms in Eqs. ((2]) and (J3]) is multiplied by 
I/ -3 . The boundary conditions are also invariant under this transformation, provided that 
the ratios among the lengths of the connectors and the perimeter of the loop are kept 
unchanged. Therefore, for this simple model we can limit our study to a single value of L, 
and the solutions for any other value are obtained by scaling. We note that L 2 ip, $ and 
Le 1//2 /£(0) = L/£(T) are invariant under this scaling. This scaling breaks down if L is not 
large in comparison to the average diffusion distance between inelastic collisions.- 14 ’ 1 " 

B. Results 

Equations (|2])-([3]) were solved numerically, as described in the Appendix. In the case 
u>s = % we studied the range L 2 /£ 2 < 10. The solutions that we found in this range 
always converged to a stationary regime. We are interested in potential differences that are 
continuous functions of <f>, but the winding number of i[} around the loop cannot change 
continuously unless -0 vanishes at some point. For this goal, sufficiently large currents are 
required. Currents that lead to appropriate behavior were found empirically. 

The lower curve in Fig. [2] shows the voltage between the electrodes as a function of $ 
for L 2 e = 10 and L 3 j = 120. The considered current density is about twice the nominal 
“critical current density” 2w(e/3) 3 / 2 ~ 70L~ 3 . Due to periodicity and symmetry of 
the interval 0 < $ < 0.5<f>o would suffice to describe this function for arbitrary flux, but a 
larger interval is presented in order to show continuity. For $ = 0.5$o; th e numeric value 
of |^| at the right extreme of the loop is of the order of 3 x lO -3, 0s; this value decreases 
if a denser computational grid is used. These and the following results were obtained for 
connectors of length 0.08L. 

Let us now consider the case w$ ^ % and denote the cross section within the loop by 
w(x). The TDGL equation and Ohm’s law become^ 

dt'ip = V 2 ijj + (e — \ip\ 2 — iip) ^ + w^w'Vifj (4) 

and 

d x <~p = u Im C0*£h/O - jw s /w , (5) 

where w' = d x w and j denotes either the current density along the connector to the super- 
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FIG. 2: Potential difference between the electrodes as a function of the flux for temperature 
(1 — 10£ 2 (0)/-L 2 )T c . is the elecrochemical potential at the normal electrode. The blue line is 
for ws = »n and the red line for ws = 3 u>n- In all cases the perimeter of the loop is 4L and 
the length of each connector is 0.08 L. For ws = u>n the current density along the connectors is 
120jo£ 3 (0)/L 3 ; for ws = 3u;n, the current density along the left connector is 75jo£ 3 (0)/T 3 . For the 
purpose of comparison between the two curves, the curve for ws = 3u>n was lowered by 15£ 2 (0)<^ 0 . 

conducting electrode, or the current density at the left extreme of the appropriate branch 
of the loop. For simplicity, we take w(x ) as the linear interpolation between ws and %. 
The upper curve in Fig. [2] shows the voltage between the electrodes as a function of <f> for 
L 2 e = 10 and L 3 j = 75 (the current density at the connector to the normal electrode is 
three times larger). In order to enhance visibility, this curve has been lowered by 15 units. 
Again, no discontinuity is visible in this curve. 

In the case ws > % we have also looked into larger values of L 2 e. Figure [3] shows the 
voltage between the electrodes for L 2 e = 30, for ws = 3wn and for ws = 6w n- To enhance 
visibility, the curve for w s = 6w^ was lowered by 40 units. 

Close to <F ~ 0.14<Fo the curve for ws = 3u>n has a steep slope, but is continuous and 
reversible. By means of an appropriate bias, this steep slope can serve to attain a high flux 
sensitivity. 

The trend suggested by Figs. [2j{3]is that the flux-modulation of the voltage increases with 
L 2 e, and the slope dip^/dQ becomes more uniform when ws/w n increases. 
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FIG. 3: Potential difference as a function of the flux for temperature (1 — 30£ 2 (0)/L 2 )T c . The 
red line stands for w s = 3wn and j = 320jo^ 3 (0)/L 3 ; the black line, for ws = 6u;n and j = 
230jo^ 3 (0)/L 3 . For the purpose of comparison between the two cases, the curve for ws = 6u>n was 
lowered by 40£ 2 (0)(/?o- The length ratios are the same as in Fig. [2j The inset shows the region near 
the inflection point for ws = 3w^ in an expanded scale. 

IV. MORE REALISTIC MODELS 

In this section we will successively refine TDGL into a realistic model. 

A. The Kramer—Watts-Tobin (KWT) model 

A model that can be justified as long as there is local equilibrium is the Kramer-Watts- 
Tobin model,—— in which the difference in the relaxation times of the absolute value and 
of the phase of the order parameter is taken into account. In this case Eq. (J4J) has to be 
replaced with 

h (d t + i(p + UT ™ ^ ^ T>i> + (e - \^\ 2 ) ^ , (6) 

where T m is the inelastic collision time and 

h = (1 + ur^\ 2 )~ 1/2 . (7) 
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FIG. 4: Potential difference as a function of the flux according to the Kramer-Watts-Tobin model, 
for a circuit with perimeter 400£(0) and ws = 3tun> for surrounding temperature 0.997T C , current 
j = 3.2 x 10 -4 j'o and ur? n = 10 4 tQ. The solid red line was evaluated assuming that the current does 
not produce heating, whereas the black line assumes a local temperature described by Eq. m 
with rj = 3 x 10 3 . The dashed line is a horizontal expansion of the solid red line, in the range where 
it has a steep slope; its vertical scale is common to the other lines, but its horizontal scale is given 
in the upper axis. 

With some algebra,— this equation can be cast in a form that is appropriate for Euler 
iterations: 

dt'i/j = h~ l \D 2 + (w'/w)V]ip — iipt/j 

+hil>{e - |-0| 2 - UT 2 n Re('ijj*['D' 2 + (w' / w)V\ij))} . (8) 

In the limit r in —> 0, KWT reduces to TDGL. 

Equations (E])-® do not obey the scalings that we found for TDGL, forcing us to fix the 
length of our circuit. 

The lower curve in Fig. Q] shows the voltage as a function of the flux for parameters that 
correspond to the lower curve in Fig. El but this time the evolution of the order parameter 
followed Eq. (jHJ) with ur 2 n = 10 4 t§, which is a typical value for low-T c superconductors with 
strong coupling. We note that the main effect of KWT in comparison to TDGL is that the 
range of fluxes for which the voltage rises steeply is shifted towards $ = O.S'Lo- We remark 
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that the flow pattern converges to a stationary situation. 

At the inflection point, d|</?N|/d( < F/ < Fo) ~ 0.35, i.e., the flux sensitivity is O^S^o/J^o- For 
T c ~ 10 K, the sensitivity is ~ 4 x 10~ 4 V/$o. Assuming that the trend is the same as in 
Fig. [21 higher sensitivities would be found for smaller ratios wg/% and for smaller lengths. 

B. Heating due to current flow 

We assume that inelastic scattering lengths are short with respect to L, so that a local 
temperature, common to electrons and phonons, can be defined. Using the units in Eq. (JT]) , 
the power deposited at a given position can be written as^ 

density of power input = ^^{j^ + uKdt + iip)^] 2 + uw~ 1 d x Re[w(Vi(;)*(dt + i(f)'ip]} . (9) 

Here j'n = jws/w — nIm(^*P^) is the normal current density, so that the first term in 
Eq. ([9]) stands for the Joule heating. The second term arises from the relaxation of the 
order parameter, and the third term from the heat carried by the supercurrent. 

Assuming that a stationary situation is achieved, this power has to diffuse away. Heat 
can either diffuse along the wire or to the substrate, but assuming that the wire is thin and 
the substrate has good thermal conductivity, heat will mainly diffuse to the substrate, and 
its flow rate will be proportional to the difference between the local temperature and that of 
the substrate; therefore, as a crude model, we estimate that the local temperature increment 
will be proportional to the density of power dissipation, i.e. 

e = e 0 - rj{j £ + u\(d t + | 2 + uw~ 1 d x Re[w(Vijj)*(dt + , (10) 

where eo corresponds to the temperature of the substrate and is a constant, expected to 
decrease when the thickness of the circuit decreases and when the thermal conductivity of 
the substrate increases. 

The upper curve in Fig. []] shows the voltage as a function of the flux for the same 
parameters as for the lower curve, but now current heating was taken into account, assuming 
Eq. (gop , with rj — 3 x lO 3 )^ 2 . We see that the effect of current heating is similar to that 
of decreasing the ratio w^/w s ; the steep rise becomes gradual. The similarity between the 
two cases may be attributed to the reduction of the ability to carry current in the vicinity 
of the normal electrode. 
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As expected, the hottest region in the circuit is found in the connector at the right, close 
to the branching point. Along the branches of the loop, the temperature gradually increases 
as the right connector is approached; this temperature is larger in the branch where the 
total current is smaller. Surprisingly, higher temperatures are encountered for <f> « 0 than 
for $ « 0.5<f>o, but in the former case the hot regions along the branches of the loop are 
shorter. 

For values of 77 that are larger (respectively smaller) than 77 = 3 x lO 3 )^ 2 , we may 
expect to obtain a 7 ?n($) curve that is similar to the upper curve in Fig. [4] provided that the 
surrounding temperature is lowered (respectively raised), so as to yield the same temperature 
at the right extreme of the loop. 

C. Thermal noise 

In this section we take thermal noise into account, and investigate to what extent it 
limits the flux resolution of the proposed device. Thermal noise affects the evolutions of 
the electrochemical potential and of the order parameter. We will add its influence using a 
formalism in which length and time are discretized. 

In the case of the electrochemical potential, we add the Johnson noise. If tpk and (pk+i 
are the electrochemical potentials in two consecutive cells, with a distance i between their 
centers, at periods of time r we have to add to p>k+i ~ pk a fluctuating term with gaussian 
distribution, zero average, and variance 

([A(^ +I - Vt )f) = , (ii) 

where Y^ = 7 re 2 £/hwa, w stands for the cross section between the two cells, and A stands 
for the deviation from the value that would be obtained in the absence of fluctuations. 

The fluctuating additions to the order parameter have been discussed in previous 
studies.— If 7 / 1 k = iV’fel ex P(b\fc) is the order parameter in cell k, then at intervals of time r 
we add to |t/A| a fluctuating term with gaussian distribution, average 

<A|*|) = EVdf-Tf , (12) 

2\^ k \ T c t 0 

and variance 

((AW) 2 } - (AWI ) 2 = r*/Af . (13) 

-L c £(3 
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where Y^ = Y (0)/ ud 2 , hk is obtained from Eq. (J7|) by setting if} = and the cross 
section w has to be taken as the average in cell k. Finally, the addition to the argument is 
a fluctuating term with gaussian distribution, zero average, and variance 


<(Aw) 2 > 


Ty, T r 
h k \^k\ 2 T c t 0 ‘ 


(14) 


Obviously, the influence of thermal fluctuations becomes negligible for sufficiently low 
temperature. Therefore, a more relevant question is whether the signatures encountered in 
the previous sections are still encountered when the temperature is lowered away from T c . 
With this in mind, we have studied the influence of thermal noise for T = 0.9 T c . Accordingly, 
we reduced the perimeter of the circuit and increased the current density. 

Figure [5] shows typical temporal variations of the voltage across the circuit, for some 
particular runs. The colored curves correspond to j = 0.095. The curve for $ = 0.2254*0 is 
typical of the case 0 < 4* < 0.44* 0 ; in this case the voltage fluctuates around some average 
value. On the other hand, for 0.44>o < 4> < 0.64>o, as seen in the curve for 4> = 0.554>o, it 
is clear that fluctuations are not normally distributed. Instead, there are large fluctuations 
that persist for long periods of time, suggesting that the circuit attains metastable regimes. 
These large fluctuations can be associated with changes in the winding number of the order 
parameter around the loop; these changes are also detected in the fraction of the total 
current that flows through a given branch of the circuit. The curves for 4> = 0.44*o and 
for 4* = 0.54>o have a shorter time span, and depict special situations for which a large 
fluctuation was present during most of the sampled time. 

When the the current density is raised to j = 0.1, large fluctuations become quite frequent 
and the flux range where they are present becomes broader. The black curve in Fig. [5] shows 
that these fluctuations are not rare for 4* = 0.254* 0 . 

It might be objected that if the system wanders between two states, it is not strictly in a 
“stationary” regime. Note however that thermal noise is present also in a standard SQUID, 
in addition to its inherent oscillatory behavior. 

The red lines in Fig. [6] show average voltages as a function of the flux for j = 0.095, 
obtained by monitoring the voltage during a lapse of time during which the winding number 
changed a handful of times. In order to assess the error, the sampling time was divided into 
49 intervals, the average voltage was evaluated for each interval, and then the statistical error 
was estimated as the standard deviation of the voltage values for these intervals, divided by 
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FIG. 5: Potential difference as a function of time for miscellaneous runs. The colored curves are 
for j = 0.095jo, and the black curve, for j = 0.1 jo- The number next to each colored curve is the 
flux in units of <ho- For = 0.225<l>o, for = 0.55d>0) and for the black curve, time is shown in the 
lower horizontal axis and the binning is 2400 to j for 4> = 0.4<f>o and for = 0.54>o, time is shown 
in the upper axis and the binning is 120 to- Other parameters: T = 0.9T C , L = 20£(0), ws = 3tCN, 
UT? n = 1Q 4 1q , rj = 0, £(0) = 10 nm, aw s = 10~ 6 fl _1 cm. For visibility, the curve for = 0.2254>o 
has been lowered by 0.1</?0) those for = 0.4<3?o and = 0.54 , o have been raised by O.l^oj and the 
curve for j = O.ljo has been raised by 0.3</?o- 

7. As a rough picture of the flux dependence of the voltage, we also evaluated the voltage 
for several fluxes, for much shorter sampling times. The results are shown as blue dots in 
Fig. [6j 

The purple lines in Fig. |H] show average voltages for j = 0.1. It appears that in this 
case our samplings involved a sufficient number of large fluctuations and the results are 
statistically meaningful. 

We can estimate the flux resolution of our circuit as 0.5<hoA</3\/At/|92(0.5 < l>o) — <p(0)|, 
where Aip is the largest statistical error of ip and At is the sampling time, and the power 
dissipated by the circuit is vj$j ||■ For the upper curve in Fig. [6j the resolution is of the 
order of 20 < F 0 to //2 - 
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FIG. 6: Average of the potential difference with thermal fluctuations, as a function of the flux. 
The red lines and the blue dots refer to j = 0.095j'o; the purple lines refer to j = 0.1 jo. The lines 
were obtained from sampling times that extended during 1.176 x 10 5 to) and represent error bars. 
The dots were obtained from samplings that lasted for 5.88 x 10 3 to> and do not show error bars. 
The gray lines are guides for the eye. The other parameters are as in Fig. [5j 

V. DISCUSSION 

We have looked for a superconducting circuit that attains a stationary regime and gives 

rise to a flux dependent voltage. We have found that this situation is indeed met within 

a very broad range of parameters. We have not conducted a systematic study to obtain 

the optimal parameters, but have rather limited ourselves to show a proof of concept for a 

proposed fluxmeter. Since the working state of this flnxmeter is not periodic, it is expected 

to exert less back action on the measured system. 

The circuit considered in Fig. [6]has a perimeter of 80£(0), which for £(0) ~ 10 nm would 

be of the order of a micrometer. For flux changes of 0.5$o, the voltage changes by about 

0.4/c#T c /e. For T c ~ 10 K this voltage is of the order of 0.3 mV, which should be readily 

1/2 

measurable. For T = 0.9 T c , the flux resolution of this circuit is rs j 20 $ 0 *o , which for 
T c ~ 10 K is of the order of 10~ 5< ho/Hz 1 / 2 . Taking awg ~ 10 _6 fl _1 cm, the power dissipated 
would be of the order of 10”'W. 

We have considered temperatures close to T c just for numerical convenience and in order 


14 





to use theoretically simple models. However, from the practical point of view, we expect that 
low temperatures will offer a better performance: lower thermal noise and reduced danger 
of uncontrolled hot spots. With all other parameters taken equal, comparison of Figs. [2] and 
El suggests that lower temperatures would yield larger signals and steeper inflection points. 
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Appendix A: Numerical solution of the dynamic equations 

The basic technique is straightforward: in each connector or branch of the loop, <p(;x) was 
obtained by trapezoidal integration of Eq. ([3]) or ([5]), and the evolution of -0 was followed by 
Euler iterations. In the following we describe the operations that were not entirely standard. 

1. Evaluation of P p 0 (p = 1,2) 

a. Case with no fluctuations 

Each branch of the loop was discretized into a grid of 2 7 vertices. A fast Fourier trans¬ 
formation was performed on 0( x ), and the Fourier transformed function was multiplied by 
(■ ik =F A) p , where k is the reciprocal variable of x. We then transformed back to x-space and 
obtained V p ip. 

Fourier transformation implicitly assumes that 0( x ) is periodic, and this feature is not 
obeyed for a single branch. In order to obtain a continuous periodic function with continuous 
first derivative, we proceeded as follows, (i) Only the T — 4 vertices in the middle of the 
grid had physical significance and represented the function -0( x ) in the physical domain 
—L < x < L] the two first and the two last values were used to interpolate between d x ijj at 
x = — L and x = L, and ensure its continuity in the extended grid, (ii) A linear function was 
subtracted from -0, so that we obtain a new function 0 / (rc) that has the same value at the 
beginning and the end of the grid; we then evaluated T )P 0 / and added to it D p (0 —0'), which 
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is easily found analytically. We could have obtained periodicity automatically by dealing 
with the order parameter in the entire loop, but this is not advantageous, because in this 
case d x fl> is not continuous at the joints between the branches. 

b. Case with fluctuations 

In this case is not a smooth function, and the Fourier method is not efficient. 
was evaluated by means of finite differences. 

2. Continuity at the contacts between the ring and the connectors 

At these points, three branches meet. Charge conservation implies that the total current 
flowing out from these points has to vanish. Moreover, for a dense computational grid, the 
amount of normal to supercurrent conversion in the region involved becomes negligible, and 
we require that the total supercurrent flowing out from the contact point has to vanish, 
i.e., lm('0o E^=i Th/E) = 0, where ipo is the order parameter at the contact point, is the 
order parameter along branch number n, and the operator V has to be interpreted as going 
outwards from the contact point. Assuming that the material in contact with the circuit 
is insulating, this condition generalizes tc3^f Yln =i = En=i — Wo Yln=i ^ n = 0, 
where A n is the component of the magnetic potential that goes out from the contact point, 
in the direction of branch n. 

In our situation Y^ n =\ = so we are left with Yln=i Ed/E = 0. This condition is 

discretized as Y^ri,=i (En,next — Eo)/ Ax n = 0, where next is the order parameter in the vertex 
along branch n next to the contact point, and Ax n is the distance to this vertex. From here 
-0o is obtained as a weighted average: -0o = (En=i ^wnext/ Ax n ) / Yln=i Clearly, this 

procedure ensures continuity of the order parameter at the contact points. 

Continuity of ip at the contact points is imposed by adjusting the division of the current 
between the two branches of the loop. 


* Electronic address: jorge.berger@braude.ac.il 
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